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ABSTRACT 

Recently, several studies have shown that young, open clusters are characterised by a consid- 
erable over-abundance in their barium content. In particular, D'Orazi et al. (2009) reported 
that in some younger clusters [Ba/Fe] can reach values as high as -0.6 dex. The work also 
identified the presence of an anti-correlation between [Ba/Fe] and cluster age. For clusters in 
the age range ~4.5 Gyr-500 Myr, this is best explained by assuming a higher contribution 
from low-mass asymptotic giant branch stars to the Galactic chemical enrichment. 

The purpose of this work is to investigate the ubiquity of the barium over-abundance 
in young stellar clusters. We analysed high-resolution spectroscopic data, focusing on the 
i-process elemental abundance for three nearby young associations, i.e. AB Doradus, Carina- 
Near, and Ursa Major The clusters have been chosen such that their age spread would com- 
plement the D'Orazi et al. (2009) study. 

We find that while the i-process elements Y, Zr, La, and Ce exhibit solar ratios in all three 
associations, Ba is over-abundant by ~0.2 dex. Current theoretical models can not reproduce 
this abundance pattern, thus we investigate whether this unusually large Ba content might be 
related to chromospheric effects. Although no correlation between [Ba/Fe] and several activity 
indicators seems to be present, we conclude that different effects could be at work which may 
(directly or indirectly) be related to the presence of hot stellar chromospheres. 

Key words: Stars: abundances-Galaxy: open clusters and associations: individual: AB Do- 
radus, Carina-Near, Ursa Major 



1 INTRODUCTION 

Elements heavier than iron cannot be produced by fusion reac- 
tions. Since the binding energy per nucleon decreases past the iron- 
group nuclei and because the Coulomb barrier increases, nuclei 
heavier than Fe are essentially synthesised by neutron («) capture 
reactions. Two distinct processes are responsible for these heav- 
ier nuclei; the slow (s) neutron capture process during stellar He 
burning, and the rapid (r) neutron capture proces fl the latter be- 
ing p r esumably related to Supemovae explosions (iRosswog et al.l 
ll999l : IWanaioetal.ll2002h . The j-process nucleosynthesis, which 
accounts for about half of the elemental abundances between Fe 
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' Neutron -capture reactions on iron seed can be slow and rapid relative to 
the average /3-decay timescale of the radioactive nuchdes. 



and Bi (see IBusso et alj[l999l : iKappeler et alj|201 ll for extensive 
reviews), is in turn produced through two channels: 

(i) the weak component, which takes place in massive stars 
(>8 Mq) during their a dvanced evo lutionar y stages of He-core or 
C-shell burning (e.g.. iPrantzos et al.lll990l : iRaiteri et al.ill993h . is 
responsible in the Solar System for the p roduction of the light est 
isotopes (atomic mass number A<90; e.g.. lTravaglio et aj]|2004) . 

(ii) The main i-process component is ascribed to thermally 
pulsing low-mass (M< 4Mq) asymptotic giant branch (AGB) stars 
and accounts for the h eaviest i-process elements, from Ba to Bi 
(e.g.. lBussoetal.ll 19991) . 

Since the abundance pattern of each component is depen- 
dent on a complex web of nuclear reactions, neutron fluxes, stel- 
lar masses and temperature ranges, i-process elements have long 
since been recognised as key diagnostics of the internal structure 
and the mixing phenomena in AGB stars. Furthermore, they pro- 
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vide excellent tracers of the chemical enrichment mechanisms in 
the Galaxy and beyond. For this reason, a wealth of observational 
studies have been performed over the years to derive i-process ele- 
ment abundances in bo t h field stars (halo, thick and thin disk, e.g 



Burris et al 



McWilliam eT^ll995l : iMashonkina & Gehrenll200l| - 
200f?) and stellar clusters (globular and open clusters. |James et al . 
2004; Yongetal.200 8; Mikolaitis et al.l2010l : |Pa"ncino et alj2010h . 

D'Orazi et al. (2009, hereafter D09) derived Ba abundances 
for a sample of 20 open clusters (OCs) spanning a wide range in 
age (~30 Myr - 8 Gyr), Galactocentric distance (7-22 kpc), and 
metallicity (-0.3<[Fe/H]<+0.4 dex). They found that the [Ba/Fe] 
ratio increases from a roughly solar value ([Ba/Fe]~0 dex) for clus- 
ters with age > 4 Gyr, up to ~0. 2-0.3 dex for ~500 Myr clusters . 
The anti-correlation between the Ba content and the cluster age can 
be reproduced by a Galactic chemical evolution model only assum- 
ing a higher Ba yield from low-mass AGB stars (i.e. 1-1.5 Mq) 
than that previously predicted (see D09 for further details). The 
same result was also confirmed by Maiorca et al. (20 11, Mil) who 
analysed both light i-process {Is; Y and Zr) and heavy .s-process 
elements (hs; La, Ce), concluding that the Ba enhancement is ac- 
companied by similar trends from the other ^-process elements. 

Remarkably, when dealing with OCs younger than 500 Myr, 
D09 detected a further increase in [Ba/Fe]: the three younger OCs 
of their sample, namely NGC 2516 (-110 Myr), IC 2391 (-50 
Myr), and IC 2602 (-35 Myr) exhibit a [Ba/Fe] ratio of 0.41±0.04, 
0.68±0.07, and 0.64±0.07 dex, respectively. While a chemical evo- 
lution model with enhanced Ba production can account for the ob- 
served raising trend up to -500 Myr, it dramatically fails in repro- 
ducing the young stellar clusters. D09 argued that indeed a process 
creating Ba in the last hundreds Myr of evolution is quite unlikely 
(unless local enrichment is invoked). In a recent paper, Desidera et 
al. (201 1) presented a thorough investigation of several fundamen- 
tal properties of the young solar-type star HD 61005, a probable 
member of the Argus association which also includes the IC 2391 
OC. The authors confirmed the previous findings by D09, obtain- 
ing that this star is characterised by extremely high Ba abundance, 
with [Ba/Fe]=0.63±0.06 dex. The origin of this unusual overabun- 
dance in Ba is not clear (see D09, Mil, and Desidera et al. 201 1) 
and several questions naturally arise from the observed pattern: is 
the enhancement in the Ba content shared by the majority (total- 
ity?) of young clusters populating the solar neighbourhood? Are the 
nearby young clusters characterised by a unique [Ba/Fe] value? Do 
they show any intrinsic internal dispersion? Do the other .s-process 
elements follow the enhancement in Ba? 

In an attempt to answer these questions, we present the first 
i-process element (both first-peak Y-Zr, and second-peak Ba-La- 
Ce) abundance determination of three young nearby associations: 
the AB Doradus, Carina-Nea r, and Ursa Major moving groups 
(see IZuckerman & Song|l2004l for a review). These three associa- 
tions cover an age range which is particularly critical for our in- 
vestigation, spanning fro m -500 Myr (UM a; King et al. 2003), to 
-200 Myr (Carina-Near; Zuc kerman et al.ll2()06l), up to -50-110 
Myr (AB Doradus; iZuckerman et all |2004 iLuhman et all |2005|; 
iMessina et al.ll201l1) . Moreover, the well-known star i Horologii, 
a probable evaporated member of the Hyades (see Vauclair et al. 
2008), is also included in our analysis. In a forthcoming paper 
(Biazzo et al. 2012, in preparation) we report iron-peak, a-, and 
odd/even-Z elemental abundances for the same sample stars, dis- 
cussing on more general grounds the chemical abundance pattern 
of these young stellar aggregates. 

The paper is organised as follows: in Section |2] we provide 
information on our target stars, while the abundance analysis tech- 



nique is described in Section |3] Our results are given in Section |4] 
and discussed in detail in Section |5] A summary of our findings 
ends the paper (Section|6f. 



2 SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

The sample was selected from the list of members of AB Dor, 
Carina-Near and Ursa Major moving groups, considering only GK 
dwarfs with projected rotation velocity vsin/ < 15 km s"' and 
no evidence of spectroscopic binarity. High S/N FEROS spectra 
(R=48000) are available for 11 stars fulfilling these criteria. We 
further added to the sample t Hor, a planet host and probable evapo- 
rated member of the Hyades, for which several abundance analyses 
based on high-dispersion spectra are available. The spectra of five 
stars were acquired as part of the program for spectroscopic char- 
acterization of the targets for the SPHERE GTO survey (Mouillet 
et al. 2010), the spectra of y Lep A and B are from Desidera et al. 
(2006), while the remaining ones were retrieved from ESO archive. 

Data reduction was performed using a modified version of 
the FEROS-DRS pipeline (running under the ESO-MIDAS context 
FEROS), yielding the wavelength-calibrated, merged, normalised 
spectrum through the following reduction steps: bias subtraction 
and bad-column masking; definition of the echelle orders on flat- 
field frames; subtraction of the background diffuse light; order ex- 
traction; order-by-order flat-fielding; determination of wavelength- 
dispersion solution by means of ThAr calibration lamp exposures; 
order-by-order normalisation, re-binning to a linear wavelength- 
scale with barycentric correction; and final merging of the echelle 
orders. 



3 DATA ANALYSIS 

3.1 Abundance measurements 

We determined Y, Zr, Ba, La, and Ce abundances for the Sun and 
for our sample stars. Spectral synthesis was applied to all the fea- 
tures under consideration, including isotopic shift and/or hyperfine 
structure (hfs) data as needed. We employed the MOOG code (Sne- 
den 1973; 2010 version) and the Kurucz (1993) set of model atmo- 
spheres, with the overshooting option switched off. Our analysis 
relies on a few lines that are proven to be reliable for this kind of 
study, namely the ^4398 A for Yn, the ^4209 A for Zr n, the ^5853 
A for Ba n, the /14087A for La n, and the ^4073 A for Ce n. We 
carefully checked in our sample spectra for other suitable lines, but 
our search was unfruitful because they are weak, or blended, or 
both (e.g., the Y n lines at A50S7A , ^5200A , ^5205A the Z r ii 
ones at /14161 A and /14443 A -see also lMashonkina et alj2007l . the 
^4322A and ^4333A of La n, the A4350 A and ^4562 A of Ce ii). 
Due to the relatively high-metallicity of our stars, we discarded the 
Ba II line at /^6141 A , which is known to be ble nded with a stronj 
iron line. Moreover, as discus sed bv Mashonkina~& GehrerJ ( 12001 
and IMashonkina et all ( |2007|) . the Ba ii Hne at /16496 A can be af- 
fected by NLTE effects that can be up to -0.2 dex (on average) 
for stars in the range -1.0 <[Fe/H] <+0.1 dex. We focused, hence, 
on the /15853 A feature, whi ch is quite insensitive to NLTE effects 
( IMashonkina & GehrenllioOQ) . as well as strong and isolated. 



There is no information in the literature on the hfs and/or iso- 
topic shifts of Y II, Zr ii, and Ce ii lines. However, the main isotope 
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Table 1. Stellai' parameters, metallicity, and i-process element abundances of our sample stars 



Star 


Teff logg f [Feii/H] 


[Yii/Fe] 


[Zrii/Fe] 


[Baii/Fe] 


[Lan/Fe] 


[Ceii/Fe] 




(K) (kms-l) 













AB Doradus 



HIPl 14530 


5600 


4.6 


1.6 


0.11 


0.00±0.10 


0.05±0.10 


0.20±0.15 


0.00±0.10 


0.03±0.10 


HIP82688 


6100 


4.6 


1.8 


0.15 






0.27±0.20 






TYC 5901-1109 


6200 


4.6 


1.5 


0.14 


0.04±0.12 


0.10±0.12 


0.40±0.15 


0.00 ±0.07 


0.06±0.12 


TYC 9493-838 


5450 


4.6 


1.8 


0.08 


0.04±0.10 


0.10±0.09 


0.30±0.15 


0.00 ±0.07 


0.06±0.14 


TYC 5155-1500 


5800 


4.6 


1.9 


0.11 


-0.05±0.10 


0.04±0.10 


0.10±0.15 


0.00 ±0.10 


0.00±0.10 


Average 










0.00±0.04 


0.07±0.03 


0.25±0.11 


0.00±0.01 


0.04±0.03 












Carina Near 










HIP58240 


5900 


4.6 


1.6 


0.02 


0.00±0.10 


0.07±0.08 


0.25±0.15 


0.10±0.10 


0.03±0.10 


HIP58241 


5750 


4.5 


1.7 


-0.01 


0.05±0.10 


0.10±0.10 


0.23±0.18 


0.00±0.07 


0.04±0.12 


HIP37918 


5550 


4.7 


1.7 


0.09 


-0.05±0.10 


0.00±0.10 


0.20±0.15 


-0.05±0.10 


0.00±0.10 


HIP37923 


5450 


4.6 


1.5 


0.10 


0.00±0.10 


0.12±0.10 


0.10±0.20 


0.03±0.07 


0.04±0.12 


Average 










0.00±0.05 


0.07±0.05 


0.20±0.07 


0.02±0.06 


0.03±0.02 












Ursa Major 










•yLepA 


6350 


4.3 


1.4 


0.05 


0.00±0.10 


0.05±0.10 


0.25±0.20 


0.04±0.10 


0.08±0.10 


yLepB 


5100 


4.6 


1.5 


-0.01 


0.00±0.15 


0.04±0.10 


0.15±0.15 


0.02±0.10 




Average 










0.00±0.01 


0.05±0.01 


0.20±0.07 


0.03±0.01 


0.08±0.10 












Hyades (?) 










(Hor 


6200 


4.5 


1.5 


0.15 


-0.05±0.10 


0.06±0.12 


0.14±0.15 


-0.03±0.10 


0.00±0.12 



of yttrium is Y, which has an odd mass number, but the level spht- 
ting is negligible because of the small spin and magnetic momen- 
tum of the yttrium nucleus. All four naturally occurring Ce isotopes 
are even isotopes (nuclear spin 1=0), i.e. '^^Cej'^^'Ce, ''^''Ce, and 
'*^Ce; due to the lack of hfs splitting and/or isotopic shifts, Ce ii 
lines are usually narrow. Of the five Zr isotopes, only one ("Zr) 
has an odd mass number and accou nts for about of th e Zr 

solar content; however, according to iMashonkina et al] ( |2007|) the 
hfs effects are negligible. Thus we adopted a single-line treatment 
to co mpute the abundances o f Y, Zr, and Ce, takin g the loggf values 
from'Hannafor d et al.l ( IT982h , lLiung et alj ( |2006|) . and lLawler eTal] 
( 12009) , respectively. 

Lanthanum, the first of the rare earth elements, has only one 
stable isotope ('^^'La), and owing to a non-zero nuclear spin (1=7/2), 
it is heavily affected by the hfs spl itting. The logg/ val ues for all its 
components were retrieved from iLawler et al. (1200 lb . Finally, al- 
though the Ba line at /i5853A does not experience severe hyperfine 
structure and isotopic shifts, we decided to employ the hfs data pro- 
vided bv lMcWilliamI ( Il998 ^ to gather robust result^ We adopted 
the isotopic solar mixtures of 81% for (''"'Ba +'^''Ba -l-''*'Ba) and 
19% for ('^'Ba -l-'^^Ba). Note that since this line shows negligible 
hyperfine structure and modest isotopic shifts, the results do not 
depend on the abundance ratio for even and odd isotopes. 



^ The hfs ti'eatment has the effect to de-saturate the spectral lines and, 
hence, can alter the inferred abundance from barely saturated lines 



Stellar parameters are spectroscopically optimised. T^ff values 
were obtained by zeroing the slope of iron abundances from Fe i 
lines with the excitation potential (typically ~ 100 Fe i lines for 
each star, covering a wide wavelength range from 4830 A to 7960 
A -see D'Orazi et al. 2009, 201 1 for details on the line list atomic 
parameters). Final microturbulence values (f) were obtained by re- 
moving the trend of iron abundances with the equivalent widths 
(EWs), while gravities were derived from the imposition of the ion- 
isation equilibrium condition, i.e. minimising the difference in iron 
abundances computed from Fe i and Fe ii lines (< 0.05 dex). The 
adopted atmospheric parameters are listed in Cols. 2,3,4 of Table[TJ 
while in Col. 5 we list the [Fe ii/H] abundances used as input metal- 
licity in the model atmosphere because all our scrutinised lines are 
at the first stage of ionisation. We note, however, that the agreement 
between [Fe i/H] and [Fe ii/H] is always within 0.02 dex for all our 
sample stars. 

In order to minimise the impact of the uncertainties on the 
atomic parameters, we carried out a strictly differential analysis 
with respect to the Sun, optimising the initial line list on a solar 
spectrum acquired with the same instrument (at the same resolu- 
tion) used for our sample stars. Adopting for the Sun Teffg=5770 
K, loggo=4.44, <fo=l-l km s"', and [Fe/H]o=0, we obtained the so- 
lar abundances listed in Table j2](Col. 2). Our result s are compared 
with the solar abundances bv lOrevesse et al.l d 19961) (Col. 3), with 
the values by Asplund et al. (2009) computed using 3D hydrody- 
namical model atmosphere (Col. 4), as well as with the meteoritic 
abundances (Col. 5) from CI carbonaceous chondrites taken from 
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Table 2. Solar elemental abundances 



4 RESULTS 



Element 


This study 


Grevesse et al. 


Asplund et al. 


Meteoritic 






(1996) 


(2009) 




Y 


2.19±0.05 


2.24±0.03 


2.21±0.05 


2.17±0.04 


Zr 


2.55±0.04 


2.60±0.03 


2.58±0.04 


2.53±0.04 


Ba 


2.13±0.05 


2.13±0.05 


2.18±0.09 


2.18±0.03 


La 


1.12±0.04 


1.17±0.07 


1.10±0.04 


1.17±0.02 


Ce 


1.60±0.07 


1.58±0.09 


1.58±0.04 


1.58±0.02 



iLodders et all ( |2009|) . Within the uncertainties, our values are in 
agreement with these previous works. 

In Figure [T] we provide an example of spectral synthesis in the 
wavelength regions surrounding the Y ii, Zr ii, Ba ii. La ii, and Ce ii 
lines for one of our sample stars (i Horologii); in the right-hand 
panels we zoomed on each feature reporting the best-fit value with 
its corresponding error (see Section 13. 2t . The plot clearly shows 
that because of the crowding of several relatively strong lines, the 
continuum placement is really critical and the synthesis technique 
is needed to derive reliable abundances. 



3.2 Error estimates 

Two kinds of uncertainties are related to abundances derived from 
spectral syntheses: (;) the errors owing to the best-fit determination 
and {ii) the uncertainties due to the adopted atmospheric parame- 
ters. 

The first source of uncertainty, which includes also errors in 
the continuum placement, ranges from 0.07 to 0.14 dex for Y, Zr, 
La, and Ce. Barium has a typical uncertainty of 0.15 dex, reaching 
a value of 0.2 dex for the fastest rotators (see Table[T). 

To evaluate the impact of stellar parameters (Tej , logg, and f ), 
we varied each quantity separately (leaving the others unchanged) 
and checked the abundance sensitivity to that variation; a change 
of ±60K in Tejf, ±0.1 km s"' in £, and ±0.1 dex in logg has been 
adopted as a conservative estimate (see D'Orazi & Randich 2009; 
Biazzo et al. 2011; D'Orazi et al. 201 1). All lines considered in this 
study are formed by ions dominating the concentration of their cor- 
responding species, and are either in the ground-state either in low- 
excitation levels. This implies that, for a given atmosphere, they 
are all equally sensitiv e to effective temperature and surface grav- 
ity variations (see also iMashonkina et al.ll2007h . We found that a 
change in Tejf of ±60K results in a variation of 0.02 dex in [X/Fe] 
and a change in logg of 0.1 dex leads to an abundance variation 
of 0.03 dex. Microturbulence dominates the errors due to stellar 
parameters: Y, Zr, La, and Ce experience a change of about 0.04 
dex, while Ba, which has the strongest feature lying on the shoul- 
der of the curve of growth, displays the highest sensitivity, namely 
A[Ba/Fe]=0.07 dex. 

The final errors can be obtained by summing in quadrature 
uncertainties given from these two different contributions, i.e., 

o"iot= -^""flt + ^\{pi,ramy Typical uncertainties are of about ~0. 10- 
0. 12 dex for Y, Zr, La, and Ce abundances, while Ba is characterised 
by larger total errors, related to the best fit determination and to the 
uncertainties on stellar parameters; both aspects reflect the almost 
saturated behaviour of this strong spectral feature. 



Our results are reported in Table [T] where the [X/Fe] ratios are 
given for all the elements along with the coiTesponding uncer- 
tainty due to the best-fit determination (Cols. 6-10). For the star 
HIP82688, a member of the AB Doradus moving group, we could 
derive only the Ba abundance, because the relatively high rota- 
tional velocity (vsim=15 km s"', see Table [3) makes the other 
spectral features blended, hampering the abundance analysis. For 
the first-peak i-process elements, in our case Y and Zr, we mea- 
sured the following average values for the three associations: 
[Y/Fe]=0.00±0.04 (rms) dex and [Zr/Fe]=0.07±0.03 dex for AB 
Doradus, [Y/Fe]=0.00±0.05 dex and [Zr/Fe]=0.07±0.05 dex for 
Carina-Near, [Y/Fe]=0.00±0.01 dex and [Zr/Fe]=0.05±0.01 dex 
for Ursa Major We found that the three young associations share 
a homogeneous solar composition as far as the light ^-process ele- 
ments (Y Zr) are concerned (the rms for all mean values is signifi- 
cantly lower than the errors). 

Focusing on the second-peak elements La and Ce, we ob- 
tained: [La/Fe]=0.00±0.01 dex and [Ce/Fe]=0.04±0.03 dex for 
AB Dor, [La/Fe]=0.02±0.06 dex and [Ce/Fe]=0.03±0.02 dex for 
Carina-Near, and [La/Fe]=0.03±0.01 dex and [Ce/Fe]=0.08±0.10 
dex for Ursa Majo]0. The hs process elements La and Ce as 
well as the first-peak i-process elements Y and Zr, show solar 
ratios. Similarly, the field star l Horologii exhibits solar abun- 
dances, with [Y/Fe]=-0.05±0.10 dex, [Zr/Fe]=0.06±0.12 dex, 
[La/Fe]=-0.03±0.10 dex, and [Ce/Fe]=0.00±0.12 dex (see Ta- 
ble[3. 

On the other hand, the mean barium abundances are 
[Ba/Fe]=0.25±0.11 dex for AB Dor, [Ba/Fe]=0.20±0.07 dex for 
Carina-Near, and [Ba/Fe]=0.20±0.07 dex for Ursa Major. This 
means that the three clusters show an enhancement in the [Ba/Fe] 
ratio of ~0.2 dex, which is not revealed in the other i-process ele- 
ments, neither the first-peak nor second-peak. In Figures |2]and|3]we 
show the [X/Fe] ratios for all the analysed elements as a function 
of effective temperature (Teff) and metallicity ([Fe/H]): with the ex- 
ception of Ba, all the i-process element abundances are consistent, 
within the uncertainties, with a solar pattern. 

The [Ba/Fe] seems also to display the largest variation within 
each association; however, since the scatter is still within the ob- 
servational errors, we argue that those variations in the measured 
[Ba/Fe] ratios are simply due to uncertainties in stellar parameters 
that are dominated by mic roturbulence. 

For the star y Lep A. lEdvardsson et al.1 ( [l993 l) -hereafter E93- 
determined Y, Zr, and Ba abundances, finding [Y/Fe]=-0.06 dex, 
[Zr/Fe]=-0.03 dex, and [Ba/Fe]=0.05 dex; their results are in very 
good agreement with our estimates for Y and Zr, while for Ba we 
obtain higher abundance (A([Ba/Fe])=0.2 dex). Differences in ef- 
fective temperature and gravity are negligible, being Teff=6398K 
(A(our-E93)=-48K) and logg=4.29 (A(our-E93)=0.01 dex); on 
the other hand they derived a notably higher microturbulence value, 
i.e. ^=1.83 km s"'. The difference of 0.4 km s"' between our study 
and that of E93 accounts for the difference in the [Ba/Fe] ratio; if 
we adopt the E93 value of ^=1.83 km s"' then we reproduce their 
derived value of [Ba/Fe] =0.05 dex. Concerning [Y/Fe] and [Zr/Fe] 
ratios, we infer values of -0.13 and -0.10 respectively. These are 
consistent with the findings of E93 (within the uncertainties). We 

' Due to the intrinsic weakness of the Ce line, we could not derive the Ce 
abundances for y Lep B (the coldest star of our sample) . We indicate as an 
average value for the cluster the [Ce/Fe] ratio measured in y Lep A quoting, 
as average eiTor, the uncertainty on the best-fit determination. 
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Figure 1. Example of spectral synthesis of Y ii, Zr ii, Ba ii. La ii. and Ce ii for the star ( Horologii. In the left panels we show the wavelength regions around 
the .v-process element features, while in right-hand panels we show the best-fit with the corresponding uncertainty. 



point out that while the microturbulence value from E93 was de- 
rived from a relationship between ^, logg and Tej (obtained from 
the analysis of about 12-17 Fe i lines only), our values are spectro- 
scopically optimised on more than 100 Fe i hues, covering a wide 
range in EP, wavelength and strength. It should be noted that in 
analysing another Ursa Major member (HR 2047), E93 derived a 
[Ba/Fe]=+0.25 dex and concluded that their Ba abundance deter- 
mination might be considered normal for young stars in the light of 
an inverse [Ba/Fe] correlation with age. Similar conclusions were 
reached bv lCastro et al.l ( [r999l) . who determined Ba abundances of 
seven GK members of the Ursa Major moving group; none of their 
stars are included in our sample, so that a direct comparison of 
stellar parameters and abundances can not be performed. However, 
their results are consistent with our findings, with [Ba/Fe] ratios 
exhibiting an excess at about ~0.3 dex level. 

The enrichment of only Barium in relation to the solar s- 



process pattern is further highlighted in Figure^] Here, we compare 
the spectrum of the solar analogue (Tefi'=5750K, logg=4.5, and 
[Fe/H]=-0.01 dex), with that of the Sun. The latter is convolved to 
the projected rotational velocity of the former. Along with Ba, we 
focused on the Zr ii and La ii, taking them as representatives of Is 
and hs elements, respectively. 

For HIP58241, we obtained a microturbulence value of ^=1.7 
km s"' and hence a difference of -1-0.6 km s"' compared to the value 
derived for the sun (£q= 1 . 1 km s"' , see Section|3]l. It can be shown 
trivially that such a difference in £ can not account for the difference 
in EWs of the Ba n line in the two stars. Given that EW(Ba n)Q = 63 
mA and EW(Ba ii). = 86 mA then (A(EW)=23 mA). If we approx- 
imate the spectral line as a rectangle (with side l=EW/2) and keep 
all the parameters constant, then we can assume that the £ values 
are proportional to the EW. Thus, a changing of +0.6 km s"' in £, 
results in an enhancement in the EW of about 12 mA. This implies 
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Table 3. Rotational velocities (vsini) from the literature and from our anal- 
ysis. 



Star 


\ sini 


\ sini 


References* 




(Literature) 


(This work) 






km s^^ 


km 




HIPl 14530 


6.6±1.2 


6.0±2 


[1] 


HIP82688 


15.0±1.5 


15.0±2 


[2] 


TYC 5901-1109 


6.0±2.0 


5.0±2 


[3] 


TYC 9493-838 


3.0±0.1 


3.0±2 


[1] 


TYC 5155-1500 


9.0±2.0 


8.0±2 


[3] 


HIPS 8240 


5.2±1.2 


5.0±2 


[1] 


HIP58241 


9.0±1.2 


9.0±2 


[1] 


HIP37918 


6.3±1.2 


6.0±2 


[1] 


HIP37923 


3.2±1.2 


3.0±2 


[1] 


yLepA 


7.7±1.8 


8.0±2 


[4] 


yLepB 


2.8±1.8 


3.0±2 


[4] 


L Hor 


6.2±1.2 


6.0±2 


[I] 



)2009l) : 14l lDesidera eTal] )2O06l) 



[3] Ida Silva et all 



that in order to reproduce the difference in the EW strength only 
in terms of microturbulence values, we should assume a value of 
^=2.4 km s"' (A(^ - ^o)=-l-1.3 km s"'). Indeed a spectral synthe- 
sis with the latter value of microturbulence yields a [Ba/Fe]=0.00. 
However, such a high value would naturally lead to a significant 
under-estimate of the iron content (with the condition of no trend 
between abundances from Fe i lines and EW strength no longer 
satisfied); similarly, abundances for the other j-process elements 
would be shifted towards lower values of about ~0. 15 dex. 

The comparison shown in Figure |4] clearly demonstrates that 
there is no difference in the strength of the Zr and La features (and 
hence in the abundances), while the Ba line is significantly stronger 
in the star. This means that the high Ba abundance we found for 
our sample stars is not affected by major systematic uncertainties, 
reinforcing the reliability of our analysis. The implications of such 
a peculiar trend in the Ba content are discussed in Section|5] 

As a by-product of the spectral synthesis technique, we de- 
rived rotational velocities (v.smi) of our sample stars. The values 
we obtain are in very good agreement, within the uncertainties, with 
previous determinations (see Table |3). 



5 DISCUSSION 

In this section we discuss the scientific implications of the derived 
abundances. Our observations suggest that the three young associ- 
ations have a solar .^-process abundance pattern concerning Y, Zr, 
La, and Ce, while Ba seems to be considerably over-abundant. In 
Section 15.11 we discuss the global .s-process element pattern ob- 
served in young open clusters and associations, while because of 
its outstanding nature, Ba is discussed separately in Section [5!2l 



■a 




40B6 4086 4067 

Wavelength (A) 

Figure 4. Comparison between the spectrum of the star HIP58241 (solid 
line) and the solar one (dotted line) convolved with the rotation profile at 
the target vsini. 



5.1 The .s-process content of young stellar clusters 



Recently, Mil derived abundances of Y, Zr, La, and Ce for 19 
OCs which complemented the Ba determinations of D09. Mil sug- 
gested that younger clusters are characterised by a higher .s-process 
element content (both first-peak and second-peak .s-process ele- 
ments) than older ones. The authors claimed that a steep increase 
is seen in OCs with ages between 1.5 and 0.5 Gyr, after which the 
[X/Fe] values remain constant around ~0.2 dex. They concluded 
that such a trend is best explained if low-mass AGB stars (1-1.5 
Mq) are more effective in the production of i-process elements than 
previously thought (see Mil for details on the anti-correlation be- 
tween the AGB mass and the efficiency in the n-capture nucleosyn- 
thesis). 

Moreover, Mil highlighted that at any given age a scatter in 
the heavy-element content is present, suggesting that it might be 
partly related to uncertainties in the [Fe/H] determinations, but that 
an intrinsic (real) variation cannot be excluded. Currently, a theo- 
retical explanation for such a scattered trend is not in hand. 

In Figure|5]we show the [X/Fe] ratios derived by Ml 1 (starred 
symbols) as a function of age (see Magrini et al. 2009), along with 
abundances for the three associations studied here. Firstly, we note 
that we can not gather reliable indication on trends of Zr and La 
with cluster age because, for those elements, only six clusters were 
analysed by Ml 1 ; consequently we focus our discussion only on Ce 
and Y. In fact, by considering only the Mil targets, the linear cor- 
relation coefficient r for Zr and La as a function of age is -0.16 and 
-0.29, respectively, and is not statistically meaningful. Conversely, 
[Y/Fe] and [Ce/Fe] ratios show an anti-correlation with the cluster 
age at more than 99% significance level, being ry=-0.59 (16 de- 
grees of freedom) and rcf-=0.66 (14 degree of freedom). Note that 
[Ce/Fe] measurements for the two youngest clusters of the Mil's 
sample are not available. If we omit Berkeley 29 (an outlier clus- 
ter, see the discussion in Ml 1) from our calculations, then we find 
ry=-0.82 and rcf=-0.59. We also computed the Speaman rank 
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Figure 2. [X/Fe] ratios as a function of T^jj for all our sample stars. Empty squares, triangles and circles refer to AB Doradus, Carina-Near, and Ursa Major, 
respectively, while ( Horologii is shown as a pentagonal symbol. Solid line represents the solar value. 



correlation coefficient (which is less affected from the presence of 
outliers) and found rspcar=-0.57 for Y (-0.68 without Be 29) and 
rspea,=-0.41 for Ce (-0.45 without Be29). 

The inclusion of our three young associations is of particular 
relevance in this context allowing us to enlarge the sample in the 
critical range from ~50 - 500 Myr. It is worth emphasising that de- 
spite the uncertainty affecting the age determination of these young 
moving groupfl the central point of our discussion is the relative 
difference in age. Thus, we can confidently assert that AB Doradus 
and Carina-Near are significantly younger than the Hyades OC, and 
slightly older than IC 2602/IC2391 (-35-50 Myr), while the Ursa 
Major moving group is marginally younger than the Hyades. 

If we consider our inferred abundances along with those of 
Ml 1, two possible scenarios can be deduced from Figure|5] 

(;) There is a systematic offset between our study and Mil. 
This may be due, for instance, to the different adopted techniques 



The age of AB Doradus is still under discussion, ranging from 50 to 
110 Myr; for our purpo se we adopt an average nominal age of ~ 70 Myr 
iZuckerman et al.l201ll) . 



(i.e., EWs vs spectral synthesis; see Section[3}. In this sense a re- 
analysis of their sample clusters by means of spectral synthesis 
would be necessary to verify this point (see also discussion in the 
following Section). 



(ii) Alternatively, we detected a true variation (at roughly 0.2 
dex level) in the .s-process abundances of young nearby stellar pop- 
ulations, co-existing both clusters with solar and over-solar heavy- 
element abundances. Under this assumption, our findings confirm 
that a scatter in ^-process elements exists for young, open clusters 
agreeing with the spread found in older clusters by Mil. We re- 
call again that there is no current explanation for such a spread (see 
also Mil). However, in order to definitely disentangle the source of 
this variation we need to enlarge the sample of .s-process element 
abundances in young OCs: a homogeneous, extensive observational 
survey on clusters younger than the Hyades is essential to draw fi- 
nal conclusions and to provide a strong observational evidence and 
constrain (perhaps new) theoretical models. 
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Figure 3. [X/Fe] ratios vs metallicity. Symbols as in Figure|2] 



5.2 The Barium issue 

In Figure |6] we plot [Ba/Fe] as a function of the age for our as- 
sociations and those OCs previously analysed by DO^|f|- within the 
uncertainties (typical values are ~ 0.15 dex), the three young asso- 
ciations fit fairly well with the global trend defined by the other OCs 
of comparable ages. The slightly lower [Ba/Fe] ratios we found for 
AB Dor and Carina-Near with respect to their coeval clusters anal- 
ysed in D09 (NGC 2516 and NGC 6475) probably reflect the fact 
that the microturbulence values are on average larger for our sample 
stars. Further evidence on the peculiar nature of the Ba abundance 
comes from the Hyades candidate-member l Horologii: for this 
star we inferred a [Ba/Fe]=0.14±0.15 dex, which is only slightly 
smaller than the average value found by D09 for the Hyades OC 
(i.e. [Ba/Fe]=0.30±0.05): given the uncertainties, we can conclude 
that the two abundances are in agreement. 

An analogous value for the Ba conten t of the Hyades 
was obtained by ICarrera & Pancinol (|20T l]), who measured 



' Because a systematic offset between clusters analysed using dwarfs and 
giants emerged from the D09 sample, we focus only on clusters for which 
the analysis is based on dwarf stars (see D09). 



[Ba/Fe]=0.36±0.04 dex. Interestingly, they confirm the remark- 
able nature of the Ba abundances: all the other .s-process 
elements exhibit solar ratios, with [Y/Fe]=-0.09±0.03 and 
[La/Fe]=-0.08±0.04 (note that they found almost identical ra- 
tios for th e coeval Praesepe clu ster). Similar outcomes were pre- 
sented bv Ide Silva et al. I j2006h , who confirmed that while Ba is 
enhanced ([Ba/Fe]=0.35±0.07), the other heavy element abun- 
dances resuk in [Zr/Fe]=-0.05±0.06, [La/Fe]=-0.14±0.07, and 
[Ce/Fe]=0.06±0.07 (i.e., solar, within the errors). 

The comparison of our results with previous studies confirms 
the Ba (and only Ba) over-abundance; Ml 1 proving to be the only 
exception. They found super-solar abundances for other elements 
in the Hyades (e.g., [X/Fe]=0.12±0.04 for Y and Ce). Because we 
are comparing studies that make use of different techniques, models 
and codes, we are left with large scope for uncertainty. We do not 
draw any conclusions on the detected abundance pattern; we only 
emphasise in this context the importance of homogeneous abun- 
dance determinations, carried out by means of spectral synthesis, 
for large sample of Hyades members, including both dwarfs and 
giants. 
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Most intriguing is the fact that none of the currently available theo- 
retical models can account for an enhancement in Ba, without bear- 
ing analogous (over-) production of the other .s-process elements, 
especially La and Ce (Bussoetal. 1999; Travaglio et al. 2004). 
This has led us to suggest that the super-solar Ba abundances are 
probably related to some (different?) conspiring effects. 

Firstly, the high Ba abundance might be related to the uncer- 
tainties in stellar parameters, especially the ones pertaining to the 
microturbulence values (see Section |3]2). At the moment, there is 
no consensus in the literature whether young, active stars are char- 
acterised by higher microt urbulence values due to the pr esence of 
strong magnetic fields (see lSteenbock & Holwege3ll98ll : James et 
al. 2006; Santos et al. 2008; D'Orazi & Randich 2009; Biazzo et 
al. 2011). However, we demonstrated in Section |4] that the higher 
^ value (in comparison to the Sun) inferred for one of our sample 
stars cannot account for the difference in the EW strength. The EW 
is significantly larger in the star implying an over-solar Ba abun- 
dance. In addition, in our analysis there is no obvious trend between 
the microturbulences (ranging from 1.4 to 1.9 km s"') and the de- 
rived [Ba/Fe] ratios; the same was found by D'Orazi et al. (2009) 
for the pre-main sequence clusters IC 2602 and IC 2391, where the 



Ba enhancement reaches levels of ~0.6 dex. As a consequence, al- 
though the microturbulence may have a significant impact on the 
inferred Ba abundances, it cannot be an exhaustive explanation. 

As further investigation, we checked if the unusually high Ba 
abundance is related to the chromospheric activity level. At solar 
metallicity, the effective depth of the /15853 A Ba ii line formation 
is at logTsooo = -1-9, while formation is at logTjooo = -3.3 and 
logTsooo = -4.8 for Ba ii lines at 6141 A and 4554 A, respectively. 
Given that the line formation zone is quite deep in the atmosphere, 
we should not expect large impact from the above hot chromo- 
spheres. As information on chromospheric activity is available for 
stars in our sample we can directly verify this prediction and look 
for any correlation between Ba overabundance and (i) Call H&K 
chromospheric emission, (ii) coronal emission and (iii) rotational 
velocity. The logRnK values were measured on the same spectra 
used for the abundance analysis as in Desidera et al. (2011), while 
for the components of y Lep logRHK was taken from Desidera et al. 
2006. X-ray luminosity for the targets was obtained looking for X- 
ray counterparts of our targets in the ROSAT All Sky Surveys (Vo- 
ges et al. 1999;2000); calibration to X-ray fluxes was performed 
following Hunsch et al. (1999). Individual X-ray luminosities for 
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the components of y Lap are from ISchmitt & Liefkd (|2004 . Pro- 
jected rotational velocities are from Table 3. As expected, we found 
that there is no relationship between these activity indicators and 
the Ba over-abundance. We note also that Ba is mostly ionised at 
the temperatures of our sample stars; therefore any over-ionising 
effect is expected to be almost negligible (see Schuler et al. 2003, 
2004; D'Orazi & Randich 2009 for over-ionisation/excitation ef- 
fects on e.g., Fe, Ti, Ca in young open clusters). 

An "indirect" chromosphere-related effect could be the shape 
of the stratification in temperature as a function of the optical depth, 
namely T(t), in the model atmosphere. Due to the presence of a hot 
chromosphere, one would expect a T(r) function less steep com- 
pared to that of old stars (the outer atmosphere should be heated at 
a certain extent by the upper chromosphere levels). To verify if this 
so we modified our model atmospheres code, manually varying the 
T(t). We found that the effect is negligible in the opposite direction 
(i.e., loge(Ba) is marginally higher). This is not surprising, because 
as explained above, the Ba ii line at 5853 A foms relatively deep 
in the stellar atmosphere. As a consequence, we can discard this 
explanation. 

Finally, NLTE barium abundances were also derived using 
the fitting between observed and calculated profiles of the /15853 
A line. The NLTE profile of this line was generated using the up- 
dated version of the MULTI code (Carlsson 1986; Korotin et al. 
1999), and the barium atomic model that is described in details in 
Andrievsky et al. (2009). For a correct comparison of the NLTE 
barium line profile with the observed one, we used a combination 
of the NLTE (MULTI) and LTE synthetic spectrum code. The LTE 
synthetic spectrum was calculated for a given wavelength range that 
includes 5853 A line. For the barium line, the corresponding h- 
factors (which express the deviation from LTE level populations), 
that were calculated with MULTL are included in the LTE syn- 
thetic spectrum code, where they are used in the calculation of the 
barium line source function. The lines of other species in vicinity 
of the barium line are synthesized using the input atomic data from 



valeB After comparison of NLTE and LTE barium abundance 
derived from the A 5853 A line, we reached the conclusion that 
NLTE effect is negligible for this line (about 0.02 dex). 

We do not possess a straightforward answer to the peculiar 
trend of the Ba abundance in young clusters/associations and sev- 
eral possibilities (including indirect effects related to the chromo- 
spheric activity level) remain open. The further enhancement at 
even younger ages (<50 Myr), reaching ~0.6 dex, deserves special 
attention, invoking urgent observational surveys for the determina- 
tion of Ba (and other i-process elements) for large samples of stars 
in very young clusters (age range 10-50 Myr). We plan to attend to 
this point in the near future. 



6 SUMMARY AND CONCLUSIONS 

In this paper we presented i-process Y, Zr, Ba, La, Ce abundance 
determinations in three nearby young associations, i.e. AB Do- 
radus, Carina-Near, and Ursa Major. Our results can be summarised 
as follows: 

- The three associations show solar abundance ratios for both 
first-peak and second-peak elements, here Y, Zr, La, and Ce. 

- The [Ba/Fe] is notably enhanced in all our sample stars, with 
average value around ~0.2 dex. 

- None of the currently available theoretical models can repro- 
duce such enhancement, accounting for a Ba production without 
a simultaneous enrichment in the other .s-process elements, espe- 
cially the second-peak ones (La and Ce). 

- We investigated whether the relatively high chromospheric ac- 
tivity of these young stars might play a role in the high Ba abun- 
dances we derived. However, we do not find the evidence of any ob- 
vious correlation between the [Ba/Fe] ratio and activity indicators. 
Moreover, uncertainties due to the abundance analysis (including 
stellar parameters), possible effects of the stratification in temper- 
ature of the model atmosphere, and NLTE corrections can be ex- 
cluded as exhaustive explanations. We stress, however, that other 
possible chromosphere-related effects can be at work and, at this 
stage, we cannot yet totally discard this scenario. 

We conclude stressing that our work reinforces the need for a 
large, homogeneous investigation of i-process abundances in clus- 
ters younger than the Hyades to draw final conclusions on this issue 
and provide observational constraints to new theoretical models. 
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